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An analytical solution for light propagation in the post-post-Newtonian approxi- 
mation is given for the Schwarzschild metric in harmonic gauge augmented by PPN 
and post-linear parameters j3, 7 and e. The solutions of both Cauchy and boundary 
problem are given. The Cauchy problem is posed using the initial position of the pho- 
ton Xq = x(to) and its propagation direction cr at minus infinity: a = - lim x(t). 

t—>-—oo 

An analytical expression for the total light deflection is given. The solutions for t — to 
and cr are given in terms of boundary conditions Xq = x{to) and x = x{t). 
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I. INTRODUCTION 

The goal of this note is to derive a rigorous analytical solution for light propagation in 
the gravitational field of one spherically symmetric body in the framework of parametrized 
post-Newtonian (PPN) formalism extended by a non-linear parameter for the component 
of order c^^ in Qij. Our interest to the analytical post-post-Newtonian solution for light 
propagation is caused by the inability of the post-Newtonian solution to predict the light 
deflection with an accuracy of 1 /las for solar system objects observed close to Jupiter and 
Saturn [1]. The maximal errors of the standard post- Newtonian solution obtained by a 
comparison of that solution with numerical integrations of the geodetic equations for light 
propagation may in some cases attain 16 /las (the details of these deviations will be discussed 
in a subsequent report) . Because of this drawback of the post-Newtonian solution the current 
implementation [1, 2] of the Gaia baseline model GREM [3] uses a kind of time-consuming 
numerical inversion procedure to compute the light deflection effect for solar system bodies 
if the accuracy below ~ 20 /ias is requested. Although it does not seem to be a problem for 
the practical use of the model, a fully analytical solution for light propagation is certainly 
desirable. 

The geodetic equation for the light ray in Schwarzschild metric can in principle be inte- 
grated exactly [4]. However, such an analytical solution is given in terms of elliptic integrals 
and is not very suitable for massive calculations. Besides that, only the trajectory of the 
photon is readily available from the literature, but not the position and velocity of a photon 
as functions of time. Fortunately, in many cases of interest approximate solutions are suffi- 
cient. The standard way to solve the geodetic equation is the well-known post-Newtonian 
approximation scheme. Normally, in current practical applications of light propagation 
in the relativistic models the first post-Newtonian solution is used. Post-post-Newtonian 
effects have been also sometimes considered [5, 6], but in a way which cannot be called self- 
consistent since no rigorous solution in the post-post-Newtonian approximation has been 
used. Such a rigorous post-post-Newtonian analytical solution for light propagation in the 
Schwarzschild metric of general relativity has been derived by Brumberg [7, 8]. Although 
the author has considered the solution in a class of gauges introducing gauge parameters, 
this parametrization does not cover alternative theories of gravity and therefore, a post-post- 
Newtonian solution for light propagation within the PPN formalism and its extension to the 
second post-Newtonian approximation is not known. However, it is clearly advantageous 
to have such a solution within the PPN formalism. The goal of this note is to repeat the 
post-post-Newtonian solution of Brumberg [7] and to extend it for the boundary problem. 

II. DIFFERENTIAL EQUATIONS OF LIGHT PROPAGATIONS AND THEIR 

INTEGRAL 

The purpose of this Section is to derive the differential equations of light propagations 
with PPN and post-linear parameters. We also find an integral of these equations from the 
condition that the solution should represent an isotropic geodetic. 

A. Metric tensor 

It is well known that in harmonic gauge 



the components of the covariant metric tensor of the Schwarzschild solution are given by 



(3) 



1 -a 

goo — - :r- — , 
1 + a 

9oi = 0, 

g,, = {l + af6,, + ^^x'x^ (2) 

1 — a 

with 

m 

a = — , 

X 

where m — is the Schwarzschild radius of a body with mass M, G is the Newtonian grav- 
itational constant, c is the light velocity, and x = \x\ = ^ 6ij x^ x^ — \J {x^Y + (^^)^ + {x^)^ 
is the Euclidean absolute value of vector x. Expanding this metric in powers of c^^, retain- 
ing only the terms relevant for the post-post-Newtonian solution for the light propagation, 
and introducing the PPN parameters (3 and 7 [9] and the post-linear parameter e one gets 

^oo = -l + 2a-2/3a2 + 0(c-*^), 
9oi = 0, 

Qij = Sij + 2^a Sij + e (^Sij + + 0{c-^) . (4) 

In general relativity one has /9 = 7 = e = 1. Parameter e should be considered as a formal 
way to trace, in the following calculations, the terms coming from the terms in gij. No 
physical meaning of e is claimed in this report. However, this parameter is equivalent to 
parameter A of Richter & Matzner [10, 11, 12] and parameter e of Epstein & Shapiro [13]. 

The corresponding contravariant components of metric tensor can be deduced from (4) 
and are given by 

gOO^_l_2a + 2{P-2)a^ + C)(c"'), 

gij^Sij-2^a5ij+(^{A^^-e)5ij-e^^ a^ + 0{c-''). (5) 
The determinant of metric tensor reads 



^ = -1 - 2 (37 - 1) a - 2 (/5 + 2e + 67 (7 - 1)) + 0(0"^), (6) 



^=l + (37-l)a + (2/3 + 4e-l + 37(7-2))a' + C(c-^). (7) 
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Metric (4) is obviously harmonic for ^ — f3 — e — 1 since the harmonic conditions (1) take 
the form 

(V^r) =(l-7)^+((l + 7f -2,3-2£)^ + 0(c-«). (8) 



B. ChristofFel symbols 

The Christoffel symbols of second kind, defined as 

2 ^ \dxP dx- dx- J ' 
can be derived from metric (4)-(5): 

r°oo = o, (10) 

T\ = ^ + {l-P)^-^ + 0{c'% (11) 

r°a = o, (12) 

r^oo = ^-(/3 + 7)^ + 0(0, (13) 

X"^ 

r^OA: = 0, (14) 



r'fc/ = 7 (a^' hi - Sii - x^ 6ik) — 

+ (2{e-^')x'5ki-{e-2^') [x' 5u + x' S,k) - 2 e ^^^] ^ + 0{c-'). (15) 

C. Isotropic condition for the null geodetic 

From now on x" denote the coordinates of a photon, denote the spatial coordinates 
of the photon, and x = \x\ is the distance of the photon from the gravitating body that is 
situated at the origin of the used reference system. The conditions that a photon follows an 
isotropic geodetic can be formulated as an equation for the four components of the coordinate 
velocity of that photon x°': 

dx" dx^ ^ 

A being the canonical parameter, or 

2 1 

goo + -goix' + ^Qijx'x^ ^0, (17) 

where a dot denotes the derivative with respect to coordinate time t. Eq. (17) is a first 
integral of motion for the differential equation for light propagation and must be valid for 
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any point of an isotropic geodetic. Substituting the ansatz x = csfi, where fi is a unit 
coordinate direction of Hght propagation {n ■ n — 1) and s = \x\/c, into Eq. (17) one gets 
for metric (4) 



s = l-(l + 7)a+i 1^-1 + 2/3-6 + 7(2 + 87) -e(^^)j a2 + C(c-^). (18) 

This formula allows one to compute the absolute value of coordinate velocity of light in the 
chosen reference system if the position of the photon and the coordinate direction of its 
propagation /x* are given. 



D. Equations of light propagation 

The geodetic equation is given by 

d'^xf' dx" dx^ _ 

ItyF^^-^lixltx-^- ^^^^ 

Using the relation of the affine parameter A and the coordinate time t — x^ / c 

d t 1 .„r\ dx dx „ / , 

one can reparametrize the geodetic equation by coordinate time t, 

F = -c' Poo - 2 c r^o,- - rj. X' + x' (c r'^oo + 2 r%. x^ + ^ r;, x^ x'^ . (21) 



Inserting the ChristofFel symbols (10)-(15), one gets the following equations of light 
propagation in post-post-Newtonian approximation, 



X 



/ r, \ ax , . axix ■ x\ 

-(^^-^x-x) ^+2(1+7) 



+2 ((/3 + 7) + (f - e) (i • i)) ^ + 2 e 



X „ c? xix- xY' 



x" x^ 



+2 (2(1 - /?) + e - 2 f ) "'^^f + 0{c-') . (22) 



Here, for estimating the analytical order of smallness of the terms we take into account that 
\x\ — C(c). Using (18) and x • x — s'^ one can simplify (22) to get 



2 ax ^ ,^ , ax (X ■ X] 



i = -(l + 7)c' — + 2(1 + 7) 



x'^ 



+2c2 (/3-e + 27(l + 7)) -;^ + 2e 



x^ 

-.2 



+2 (2(1 - /?) + e - 2 7^) " + 0{c-') . (23) 
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E. Equations of light propagation with additional empirical parameter a 

For our purposes it is advantageous to have one more additional parameter that can be 
used to trace terms in the folowing calculations which come from the post-post-Newtonian 
terms in the equations of motion of a photon. We denote this parameter by a and introduce 
it in the above equation simply as a factor for all the post-post-Newtonian terms: 

^ r, ax , ^axiic- x\ 
^ = -(l + 7)c^^ + 2(l + 7) \^ ^ 

+2c2a(/3- 6 + 27(1 + 7)) ^ + 2ae^^^4^ 



+2 a (2(1 - /?) + 6 - 2 "'^^t'"^^ + ^^(c"') ■ (24) 

Setting q; = in the solution of Eq. (24) one can get formally a second-order solution for 
the post-Newtonian equations of light propagation. 

III. INITIAL VALUE PROBLEM 

Let us now solve analytically an initial value problem for the derived equations. 

A. Analytical post-post-Newtonian solution 

For two initial conditions 

cco = a;(to) , 

0-= hm (25) 

using the same approach as was used by Brumberg [7], one gets a formal solution x{t) and 
icit) of Eq. (24): 

-ijv^o-, (26) 
c 

ajjv = cco + c (i - to) (27) 

-cCpN = <T + m Ai(a;Ar), (28) 

x^^ = Xn + m {Bi{xf^) - Bx{xq)) , (29) 

-*ppN = o- + "^^i(a5pN) + "^^^2(a;iv), (30) 

,2 



iCppN^aJiv + m (Bi(a;pN) - Bi(a;o)) +m {B^ix^^) - B2{xq)) , (31) 



with 



Aia; =- 1+7 f + - ' 32 

\x\x — a x) X j 

Bi(a;) = -(l + 7) (^^^^-^^^ + o-log(x + o--a;)V (33) 
\ X — a ■ X I 
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1 / / s / s ^ s / . (7 X (x X cr ) 

-- 8 1+7-q;7 1 + 7 -4Q;/3 + 3Q;e <T • a; ^ 

4 ^2 |cr X a;|2 

(8(l + 7-a7)(l + 7)-4«/? + 3ae) {n-S{a,x)) (34) 



„ , , ,^ ,2 u ,^ ,0 o- X (cc X <t) 1 a; 

^2 3^ = - 1 + 7 + 1 + 7 7 ^ + 7 « e ^ 

x-crx {x-cr-x) 4 



4 a (8 (1 + 7) - 4 /? + 3 e) - <5(^t, a.; 

4 |<T X a;| Vz 

4 (8(l+7-Q;7)(l + 7)-4Q;/3 + 3Q;e) {<T ■ x) (tt - a;)) , 



(35) 



or, alternatively, for Bi and B2 



„ , , , /cr X (a; X cr) , , ,\ , , 

Bi a; 1 + 7 ^ ^-o-log x-o--a; , 36 

\ X — a • X J 

„ , , ,^ ,2 (T ,^ ,2 o- X (a; X ft) 1 X 

^2 a; = + 1 + 7^ + 1 + 7 2- ^ d^ + 7«e^ 

X - <r • a; (x - <r • a;) 4 

a (8 (1 + 7) - 4 /3 + 3 e) (| - 5((t, x; 

(8(1 + 7-q;7) (1 + 7) -4Q;/3 + 3Q;e) (<r • a;) ^ ^ (tt - 5(<r, a;)) . 

(37) 

Here 5(a, 6) is the angle between two arbitrary vectors a and h. Clearly, for an angle between 
two vectors one has < 5(a, 6) < tt. The angle (5(a, 6) can be computed in many ways, for 

example as 5 (a, 6) = arccos — — . With these definitions the solution of (24) reads 

a 

a; = a;ppN + 0{c~^), 
^a;=^a;ppN + 0(c-^). (38) 

It is easy to check that the solution for coordinate velocity of light cCppN satisfies the 
integral (18). In order to demonstrate this fact, it is important to understand that position 
X in (18) lies on the trajectory of the photon and must be therefore considered as ajpN in 
the post-Newtonian terms and as a^N in the post-post-Newtonian terms of (30). 
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B. Total light deflection 

In order to derive the total light deflection, we have to consider the limits of coordinate 
light velocity cb for i — > ±00: 

lim -x(t)^a, (39) 

c 

lim -x{t)=u 



= (T-2 1 + 7m , ^" ,„^ -2 1 + 7fm'- ^ 

1 /^/. V/- V , ^ ^ V 9 cr X (xn X cr) 

-TT (8(1 + 7-q;7)(1 + 7) - 4Q;/3 + 3Q;e) ^ ^ 



4 |a;o X cr 

X (ccq X 

\Xo X <T 



+2 (1 + 7)^ (xo + • a^o) Z^":" + ^(^"')- (40) 



Accordingly, the total light deflection reads 

|<T X 1^1 = 2 (1 + 7) m 7 -2(1 + 7)^ {xq + <t • Xq 



|3 



cco X cr| \Xq X cr| 

+7(8(1 + 7-q;7)(1 + 7) - 4:aP + 3ae) nm^ ^^ + 0(0"^). (41) 

4 |a;o X crp 

Eq. (41) defines the sine of the angle of total light defiection in post-post-Newtonian ap- 
proximation. The first term in (41) is the post-Newtonian expression of total light defiection 
~ Am/d, where d— \xqX(t\ being the impact parameter for the fight trajectory. The other 
two terms are the post-post-Newtonian corrections. Note that although the total light deflec- 
tion |cr X i/l is a coordinate-independent quantity, Xq and, therefore, the impact parameter 
d = \xo X cr| are coordinate-dependent. In order to compare Eq. (41) with results expressed 
through the coordinate-independent impact parameter d' one has to find a transformation 
between that impact parameter and d. The coordinate-independent impact parameter is 
defined as a limit 

d' — lim \x{t) X cr| = lim \x{t) x i/\ . (42) 
t— >— 00 

Substituting the post-Newtonian coordinates (29) of the photon x{t) into this definition one 
gets 

d'^d+il + ^)m ''"^J'^" + Oic-') . (43) 

It is now clear that the second term in the right-hand side of (41) just "corrects" the main 
post-Newtonian term converting it to 2(1 + ■^,)m/d'. Using d' one can write (41) as 

771 1 ?T?^ 

|<rxi/| = 2(l + 7)- + - (8(1+7-q;7)(1 + 7) - 4:aP + 3ae) tt — + 0(c-^) . (44) 

This result with a = 1 coincides with Eq. (4) of Epstein & Shapiro [13] and also agrees with 
the results of Richter & Matzner [10], Cowling [14] and Brumberg [7] in the corresponding 
limits. 
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IV. BOUNDARY PROBLEM 

For practical modeling of observations of solar system objects it is not sufficient to consider 
the initial value problem for light propagation. Two-point boundary value problem is of 
interest here. Let us consider that the light ray must propagate between two points being 
emitted at a position Xo at time moment to and received at position a; at a time moment t 

Xo^x{to) , 

x^x{t). (45) 

Initial time moment to can be considered here to be known (although for any stationary 
metric like that considered here this moment plays no role), but the final moment t is 
unknown. We also denote 

R = x-Xo, (46) 
k-^, (47) 

where i? = is the absolute value. In the following, the solution of geodetic equation (24) 
will be expressed as a function of the boundary values Xq and x. 

A. Formal expressions 

An iterative solution of (26)-(31) for the propagation time t — t — to and unit direction 
cr reads as follows: 

cT^R-mk- [Bi{x) - Bi{xo)] - k ■ [-62(0;) - B2{xo)] 

+ ^\kx {B,{x) - B,{xo))f + 0{c-'), (48) 

a-^k + m^ {kx[kx {Bi{x) - Bi{xo))]) 
R 

W - (fc X [fe X (^2(0;) - B2{Xo))]) 
K 

W {B^{x) - B,{xo)) x[kx {B,{x) - B,{xo))] 

-^m'^k\kx (B,(x) - B,(xo))f + 0(c-'). (49) 

These expressions are still implicit since for the post-post-Newtonian accuracy one should 
use the post-Newtonian relation between <t and k in Bi appearing in the post-Newtonian 
terms. That relation can be again obtained from (49) by neglecting all terms of order 
C(m^). On the contrary, in the terms in (48) and (49) of the order of C(m^) one can use 
the Newtonian relation a — k. 
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B. Propagation time 

Substituting (33) and (35) into (48) one can derive an explicit formula for the time of 
light propagation: 

X + xq + R 



CT — R + (l + 7)m log 



X -\- xq — R 



1 fxl-x^-R'' x^-xl-R^\ 

— + — i — j 

+ ] a (8(1 + 7) - 4/3 + 3e) , ^ , 6{x, Xq) 

4 |£C X CCol 

This formula allows one to compute time of light propagation r for given boundary conditions 
Xq and X. Here we have used that 6{k, Xq) — 6{k, x) = 6{x, Xq). 

C. Transformation of fe to cr 

In the same way substituting (33) and (35) into (49) one gets: 

, ,^ , X — xo + R , , , 

a = k + (l + 7)m ^ fe x (ccn x x) 

\x X Xof 

(1 + 7)^ . {x-xo + Rf 



2 \x X Xq\ 



+ m k X {xq X x) 



1 1/1 1 



JU JUr\ 



1 1 

+ - (8(1 +7 - q;7)(1 + 7) -Aa(3 + 3ae 



X X Xo 



3 



X (2R^ (tt - 5{k, x)) + (x^ -xl- R^) 5{x, Xq)^ 

+ 0{c-') . (51) 



+^ (1 + 7)2 , (x + Xo) (x - Xo - i?) (x - Xo + i?)2 



.1(1+^)2^ 1 

2^ ^' \xxxo 

This formula allows one to compute cr for given boundary conditions Xq and x 



D. Transformation of cr to n 

Considering the relativistic model of positional observations (e.g., Klioner [3]) it is clear 
that a unit direction of light propagation at the point of reception n plays here an important 
role. This vector is defined as 

x{t) . ^, 

Expanding Eq. (30) up to order 0{c ^) leads in post-post-Newtonian approximation 

n = <r + m Ci(a;pN) + C2{x^) + 0(c-^) (53) 
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with 



Ci{x)=A,{x) - a{a-Ai{x)) = -(I + 7) 



(T X {x X (t) 
X {x — (T ■ x) ' 



C2{x)^A2{x) - Ai{x) {(T-Ai{x)) - - a- iAi{x) ■ Ai{x)) - a {a ■ A^ix)) 



a ■ X 



— X (cc X 0-) + (1 + 7)' 



x^ 



+(1+7)' 



(J X [X X (J) 



X [X — cr ■ X] 



0(1 + 7)^-^ 



(7 X [X X a) 

x"^ {x — <T ■ X) 
<7 X + (T ■ X 



x^ X — a ■ X 



-- 8 1 + 7-q;7 1 + 7 -4Q;/3 + 3Q;e {a ■ x) ^ 

4 |cr X £cp 

-J (8(1 + 7-^7) (1 + 7) -4a/3 + 3ae) """^^ (tt - 5(o-, a;)) . (54) 

Noting tliat Xpj^ = x + 0(c^'^), a; = + -RA; and using Eq. (51) for a in Ci(a;pN) we get 



n — cr — (l + 7)mfex (xq x x) 



R 



X X Xq 



1 + 



k ■ X 



X 



+^(l + 7)2m2 
fc X (cco X cc) 



a; X cco r ^ 
2 



— I 1 H I (3x - xo - -R) (x - xo + R) 



X 



(i+7r 



X X Xq 

(8(l + 7-«7)(l + 7) -4a/3 +3ae; 



\ ^ fc-a;\ p (i^^ - (x - xq)^) ^ 1^ 



2\x X xo\^ 
k ■ X R 



x'^ \x X Xq\^ 



1 n- 

-- (8(l + 7-a7)(l + 7) -4a/5 +3ae) , ■ ^ (tt - 5(A;, a;)) 



+ 0(c- 



(55) 



This expression allows one to compute the difference between the vectors n and «t, starting 
from the boundary conditions Xq and x. 

E. Transformation of fc to n 



Finally, a direct relation between vectors fc and n should be derived. To this end, we 
directly combine Eqs. (51) and (55) to get 

/-. \ fc X (a^o X a;) ,^ ^. 
n = fc - 1 + 7) m — ^— '— 1 + F) 

x[xXq + X ■ Xq) 

m? {{x — xqY — R^)^ 



+ fc X (a;o X x) 



\X X Xq\ 



(1 + 7)^ 



R^ -{x- xof 



2 ^ ' " la; X a;oP 
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with 




-2 



k ■ cc' 



1 
4 
1 



(1 + 7 - q;7)(1 + 7) - 4:aP + 3ae) R 



k ■ X 



x^ 



'1 + 7 - q;7)(1 + 7) - AaP + 3ae) 



X 



X X Xq 

R^ 







\X X Xq\ 



5{x,xq) 



+ O (c-^) (56) 



-(1 + 7) m ■ 



X-\- Xq 



(57) 

XXq + X ■ Xq 

This formula allows one directly to compute the unit coordinate direction of light propagation 
n at the point of reception starting from the positions of the source Xq and observer x. 
The reason for writing the post-post-Newtonian term proportional to F together with the 
main post-Newtonian term is that it is the only post-post-Newtonian term which cannot be 
estimated as const • m^. This will be discussed in a subsequent report. 



V. CONCLUSION 

The analytical post-post-Newtonian solution for light propagation derived in this note 
will be used in a subsequent report to obtain an analytical formula for light deflection for 
solar system bodies. Not all post-post-Newtonian terms in the above formulas should be 
used to obtain the goal accuracy of I /xas. Analytical estimations of individual terms can 
be used to find those of them which are numerically relevant at the level of I /xas. Detailed 
estimations and comparison to numerical solutions of the boundary problem will be given 
elsewhere. 
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